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Abstract 



A new parameterization of unconstrained degrees of freedom for gravitational 
field, used in |S] , has been generalized to one-parameter family of such parameter- 
izations, depending on a real parameter (3 S [0,2]. The description introduced in 
|H] corresponds to the special choice (3 = 0. The method is closely related to the 
proof of the positivity of the energy presented in [3] where /3-foliations have been 
introduced (see also applications to black holes dynamics in 0], |S] and [7]). 

Spherically symmetric initial data corresponding to trivial degrees of freedom 
is analyzed along these lines. In particular, the quasi-local energy content of the 
Schwarzschild initial data is analyzed for different choices of the /3-gauge. 

1 Introduction 

We consider a compact, smooth, three-dimensional manifold V, diffeomorphic to 



Denote by dV its boundary. Limiting cases ro — ► and/or r\ — > oo may be also 
considered. 

Consider the spacetime "tube" M = V x R 1 and its boundary T = dV x M 1 which 
is a one-timelike and two-spacelike hypersurface in our spacetime. Choose coordinates 
on M in such a way that (x 1 ^ 2 ) are coordinates on dV (e.g. spherical angles 9 




and (p), x 3 = r is a ,, radial" coordinate which is constant on dV. Moreover, denote by 
x° the time coordinate. So we have 

V t := {x e M : x° = t} = [j S(r) where S(r) := {x £ V : x 3 = r} , 

r£[ro,n] 

T = {x e M : x 3 = r } U {x G M : x 3 = n} . 

We use the following convention: Greek indices fi, v, . . . label spacetime coordinates 
and run from to 3; Latin indices k,l, . . . label space coordinates on V and run from 
1 to 3; Capital indices A,B,... label coordinates on dV (,, spherical angels") and run 
from 1 to 2. 

Consider Cauchy data (gki,P kl ) for Einstein equations in the three-dimensional 
bounded volume Vt with boundary dV%. This means that g^i is a Riemannian metric on 
Vt and P kl is a symmetric tensor density which we identify with the ADM momentum 
(seed) 

P kl = s/A^Z{g k hxK-K hl ) . 

Here, K kl is the second fundamental form (external curvature) of the imbedding of Vt 

into the spacetime M which we are going to construct. 

The twelve functions (gki,P kl ) must fulfill four Gauss-Codazzi constraints: 

(det g mn )K - P kl P kl + \{P kl g k if = 167r(det g m n)T^n v , (1) 



P^l = 87T^detg mn Ti^, (2) 

where T^ u is the energy momentum tensor of the matter. By 1Z we denote the (three- 
dimensional) scalar curvature of gti , whereas is a future timelike four- vector normal 
to the hypersurface Vt- The geometric structure used in JJ) and © (the covariant 
derivative "|", rising and lowering of the indices etc.) is the one defined by the three- 
metric g H . 

Einstein equations and the definition of the metric connection imply the first order 
(in time) differential equations for g^i and P kl (see or j2] p. 525) and contain the 
lapse function iV and the shift vector N k as free parameters, canonically conjugate to 
the four constraints (JTJ) and @: 



m = ^ {p k i - \g k iP^j + N w + N llk , (3) 

where g := detg mn and P := P kl gn, 

P kl = -N^g (n kl - \g kl T^j - ^ (p km P m l ~ ^PP kl ^j + {p kl N m ^j | + 



2 



+^9 kl (p kl Pki ~ \P^J ~ N k lm P ml - N l lm P mk + (n^ 1 - g kl N\ m \ m ) + 



+MN^T mn g* m g ln . (4) 

2 Reduced phase space of Cauchy data 

We want to analyze Cauchy data in terms of the 2+1 decomposition of the initial 
surface. To describe the independent degrees of freedom of the gravitational field we 
use the following objects (w k , Si): 

w k := \g 3k (g 33 )"' 1 , (5) 
where A := y/det gAB is the two-dimensional volume form, and 

s A :=(w 3 y l P 3 A S3 :=-( w y l (^S + P 3 A 9 ^ , (6) 

where ~g AB is the two-dimensional inverse of the two-metric g AB , whereas by S we 
denote g AB P AB . As will be seen later, the data (w k ,si) describe partially reduced 
phase space of gravitational field. 

Denote by a AB the standard metric on a unit sphere S 2 , {a AB dx A (g> dx B = d# <S> 
d6 + sin 2 6 dip <g) dip) and by a = Vdet <j A b(= sin 9) the corresponding volume element 
on the unit sphere. On each sphere S(r) we introduce the following two-metric 

VAB ■= fy\~ l g AB , (7) 

conformally equivalent to g AB . Its inverse metric is given by fi AB = a- x \~g AB . 

It is easy to check that the left-hand sides of the vector constraints Q may be 
rewritten as follows: 

I / p33 \ 

PA k \k = {Xg AC S AB ) l{c + (s A w k ), k - w k s KA + - ( -33- + PS J g 33 , A (8) 

1 „3k fP 33 w k \ (w 3 g AB s B \ * Co (w B \ 

g^ P \ k =[^^)^ k+ [-^)^- X9 SbC \^) u + 

+\S AB (Xg AB ) ,3 + ^s lW k , k + \ + /3s) ^ 33 , fc (9) 

where S AB := \~ 1 {P AB — ^S) and "||" denotes two-dimensional covariant derivative 
with respect to the two- metric fi AB . 
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To describe effectively 1 the reduced phase space, i.e. the space of classes of gauge 
equivalent pairs (gki,P kl ), one is free to impose four gauge conditions which enable us 
to pick up a single representative within each gauge-equivalence class. 

We propose the following conditions: 

p33 

+ ps = , (10) 

g 66 

= 0, (11) 

fJ-AB = °AB ■ (12) 

It is easy to check that the gauge condition (|11[) may be rewritten in the following form: 

h n 3k 2 

fe(ln 5 33 ), fc -^, (13) 



v /^33 ^33^* /'« x 3 

where k is the two-dimensional trace of the extrinsic curvature kAB of S(r) with respect 

to the three-metric g\a on Vt- 

Conditions (|10j) and describe a specific "2+2" decomposition of spacetime. The 

two-parameter family of surfaces t = x° = const., r = x 3 = const, (topological spheres) 

is defined in terms of a nonlinear system of partial differential equations 2 imposed on 

coordinates r and t. More precisely, eq. (|l(Jj) rewritten in terms of extrinsic curvature 
K ki. 

l-(3K™ 
l + (5 g 33 

leads to the following PDE for the unknown functions t and r: 

/ VH \ = l-l ( V u t \ (dr\dt)V»t - (dt\dt)V»r 
" \y(-dt|dt)/ 1 + /? "{y/FdtldtjJ (dr|dt) 2 - (dt|dt)(dr|dr) ' 1 j 

where here the notation is four-dimensional with respect to the four-metric g^, e.g. 
(dr|dt) := g^^ 'd^rd v t. Similarly, the gauge condition (fTT|) takes the following form: 



[(dr|dt) 2 - (dt|dt)(dr|dr)] 
r 2 {-dt\dtf 



(drdt) „ 
(dt dt) 



(15) 



However, for (3=1 the construction splits into separate equations, the first one gives a 
maximal three-surface t = const, and the second one corresponds to a certain spherical 
foliation of this maximal surface (conformally harmonic gauge in 0). 

1 The procedure proposed here does not cover the entire phase space but only an open neighbourhood 
of the flat initial data. 

2 They are quite close to the gauge conditions discussed in |10|. 
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Condition l|12[) corresponds to the choice of appropriate conformal coordinates on 
each two-dimensional sphere S(r). This is possible because every two-dimensional 
topological sphere is conformally equivalent to a unit sphere and eq. (|12j) describes 
precisely this equivalence. However, coordinates (x A ) are not fixed uniquely by the 
above condition but only up to the six-parameter family of conformal transformations 
of the unit sphere. We have, therefore, the residual, six-dimensional gauge freedom on 
each S(r). Using this freedom we may annihilate the six-dimensional dipole component 
of the vector w k . 

Due to gauge conditions and (j!2|) we can simplify vector constraints as follows: 



<tSa B \\B + {s A W k ), k - w"s kjA = 87TJA 



(16) 



where ja ■= \/det gm n T A ^, and 
,,m„,,k 



20 



s m w w 



w 



• k + 



B 



, A ~ cjS a b ( ^3 ) + —si = 8vrj 3 , (17) 

" \\ A 



2w 
r 



where 73 := v/ det g mn Tk u n^^r. Observe that for 0=1 the "conformal factor" g 33 

does not enter into the equation. 

Equation (|16|) is a two-dimensional internal equation on each slice S(r) separately 
and it enables one to reconstruct Sab from the reduced data (sf.,w l ). There are, 
however, additional constraints which must be fulfilled by this data. These constraints 
are visible if we contract (|16|) with an arbitrary generator ^ A of the six-dimensional 
conformal group on S(r). Such a generator fulfills the conformal Killing equation 

€a\\b + €b\\A = °AB i°\\c 1 



and, whence, 



(vs A B i A ) 



\\B 



(vS A B Z A ) 



The integral of this expression over S(r) vanishes identically. Consequently, we have 
six residual constraints imposed on the reduced data (sk,w l ) on each sphere S(r): 



S(r) 



.S A W 



8vr 



t A 3A ■ 



S(r) 



(18) 



They enable us to calculate the (six-dimensional) dipole part of s A , canonically conju- 
gate to the dipole part of w A , which was annihilated by the residual gauge condition. 

Equation (|17|) provides a relation between variables s^, which are no longer inde- 
pendent parameters. If we know sa and w s 7^ than 1)1 7JI can be viewed as an equation 
for the unknown function S3 with given sa- Similarly, the gauge condition (jll|) enables 
us to calculate w 3 once we know w A . 
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We conclude that, similarly as in the (partially) reduced canonical data (w k , si), 
fulfilling the gauge condition Qll[) and the residual gauge condition (the dipole part of 



w 



vanishes) contain all the information about the complete Cauchy data (P ,gu) 
satisfying the constraint eqs. © and the gauge conditions (fTT|) . ((T2*]) . In particular, 
to calculate the "conformal factor" g 33 we must solve the scalar constraint (Q, which 
can be rewritten as follows (see also 6 ): 



,Xw k , , 

2(^*0,* 



2d, 



w 



Xg AB 



B 



+ 



/733 V 2 



A 



PkiP 



:P 2 ) + 



'9 



k A Bk 



AB 



1 



-k 2 + Wirp 



(19) 



where p := T^ iv n^n u is a matter density, k AB is an extrinsic curvature of two-surface 
S(r) and i? is a scalar curvature of induced two-metric g AB (see (3], (U and [5]). 
Let us notice that the following identity holds: 



PkiP kl ~ ^P* = ^ ( V ) + ^^a^b + -Xq^SabScd - 



1 /P 



■33 



,,33 



-AB D 3 



-AC* \ ~DB « 



3 33 



,,33 



.33 



— (w^y-g^SASB + Xg AC Xg UB S AB S C D + 2/3(j3 + 2)(w k s k Y , 



(20) 



where the last equality is implied by the gauge condition — ==- + (5S = 0. Let us observe 

that (|2*U)) is nonnegative for /3 > 0. Moreover, due to gauge condition ((TTj) and identity 
(J20J), the scalar constraint (|19|) takes the following form: 



V 3 ^ + 



g 33 \ Xw K 



, k + XR+[X~g AB (logg 33 ), B ] , A 



(w 3 ) 2 g AB s A s B + 



33 



+ X {^^SabScd + 2/3(/3 + 2)(*/s fc ) 2 ) + ^3/3(2 - (3) ( , J 



+A (k AB k AB - \k 2 + i^(log 5 33 ), A(log 5 33 ), B + levrp) 



(21) 



It is easy to check that for (3 £ [0, 2] the right-hand side is nonnegative. This observation 
enabled us to show the positivity of ADM mass (see [3]). For this purpose we observed 
that the integral of the left hand side over the whole K(0, ro,r%) for ro — > and r\ —* 00 
gives the surface term at infinity, proportional to the ADM mass. 
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3 Symplectic structure and the complete reduction 

Let us observe the following identity: 

/ p33 \ 

-P kl dg kl = 2s k dw k + aS AB dfi AB +l — + ps) dln 5 33 , (22) 

which may be easily checked by direct computation, using definitions ©, © and Q. 
It leads to the partial reduction of the phase space if we impose gauge conditions (fTU|) 
and (fT2*)). More precisely, the last two terms drop out and we are left with 

-P kl dg kl = 2s k dw k . 

Defining 




and 

2r 2 • s k = A k , 
we get the following symplectic structure: 

-P kl dg kl = A k dD k , (23) 
together with the constraint which now reads: 

d k D k = . 

This structure is formally equivalent to the structure of classical electrodynamics, where 
D k is the electric displacement vector-density and A k is the vector potential for the 
magnetic field. Constraint ()17j) plays a role of the non-linear "Coulomb gauge condi- 
tion" in electrodynamics. The analogy is not complete, because here we have also the 
residual gauge condition (dipole part of D k vanishes), dual to the residual constraints 

Further reduction of the phase space (D k , A{) may be performed if we observe that 
only two variables among D k and two variables among A\ are independent. Equations 
(|TT|) and (|17jl. together with the scalar constraint (^Q), define a subspace in the space 
of variables (D k ,Ai)), corresponding to the two degrees of freedom of the gravitational 
field. To parameterize this space by independent variables, we may represent 

Ai = Ai + di4> , 

and impose a further gauge condition: 

A B = sbc V CD d D ip . 



7 



Inserting this Ansatz into Q23|) and integrating over the entire V, we obtain the following 
symplectic form: 

n = J (dD 3 Adi 3 + d (e BC d B D c ) Ad^ . 

It is obvious that the above four functions contain the entire information about (D k , A\). 
To reconstruct the physical data (P kl ,gki) we must solve constraint equations. In 
particular, eq. (|17|) is a three-dimensional, elliptic equation for the function (f>. 

4 Spherically symmetric initial data 

Let us assume that our initial data are spherically symmetric and the gauge condition 
is compatible with spherical foliation. This simply means that sa = wa = and, 
moreover, w 3 = r 2 a which is consistent with The energy-momentum tensor is no 

longer free, we have that ja = 0, T 3 a = and, moreover, that the traceless part of Tab 
vanishes: 

Tab - -^9ab9 CD Tcd = , 
The vector constraints (|16|). (|17j) simplify drastically: 

aSA B HB-w 3 s 3A = 0, (24) 



2w 

2p(s 3 w 3 ) , 3 + s 3 = j 3 



(25) 



Assuming that 3a3z = (spherical symmetry of j'3) we deduce from (|25|) that ^^3 = 
and from we get Sab = 0. 

The scalar constraint is also very simple: 



V3fe + 



9 



X + ra 



= 2f3((3 + 2)^-(w 3 s 3 ) 2 + ^3/3(2 - (3) (<? 33 , 3) 2 + IGnXp . (26) 

4.1 Special cases (3 = 0, |, 1 and the Schwarzschild initial data 

We would like to compare different expressions for the Schwarzschild metric, which we 
obtain for different /3-foliations. 
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4.1.1 Conformal gauge (3 = 1 



Let us start with conformally flat representation of the Schwarzschild initial data three- 
metric: 

4 

„2 , 2 A A 



ds 2 = (l + ^) (dr 2 + rV^darW) 



&AB ■ 



(27) 



Obviously r = corresponds to the second spatial end. One can easily check that 



w 



r 2 a . 



n 33 



1 + 



my 
27/ 



k 



2(r 



/^33 r(r+f) 



and 



5 := 2 



,,33 



A + ra 



4ma 



2ma for r 
4m<r for r : 



m 
2 
CO 



Let us observe that the surface integral 
1 



16vr 



Sir) 



(28) 



which we obtain when integrating the scalar constraint ()26|) over V, gives half of the 
ADM mass when calculated on the minimal surface S(r = f) and the entire ADM 
mass for S(r = oo). This means that the "energy density" — the right-hand side of 
the scalar constraint l|26|) — splits into a half which is contained outside the horizon 
and another half, hidden inside the horizon. The ADM mass seen at both space ends 
is equal to each other. 



4.1.2 Harmonic gauge = 1/2 

Let us denote by R = x 3 the solution of eq. (|lljl for j3 = 1/2, to avoid confusion with 
conformal coordinate r introduced in the previous subsection. It is easy to check that 
R = r + f . This implies the following: 



33 
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9 + E R -) X + Ra 



r + — 
4ma- 



r + 



3ma 



for 
for 



r 

r ■ 



2 

oo 



Hence, the surface integral (|28jl gives 3/4 of the ADM mass when calculated on the 
minimal surface (and, again, the entire ADM mass for S(r = oo)). 



4.1.3 Inverse mean curvature flow (3 = 

This particular case has been already analyzed in and 0. The gauge condition (|1 1|) 
or, equivalently, equation 1|13|) . corresponds to the so-called inverse mean curvature 
flow (see (§]) and its solution enables one to prove the Penrose inequality (see jllj 
and references therein). In this case equation (|llj) implies that the coordinate x 3 = 

1 H I is the usual Schwarzschild coordinate. Moreover, for the Schwarzschild 

initial data and (3 = the surface integral ()28|) is constant on each sphere S{r) and 
plays a role of a quasi-local mass. 
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5 Conclusions 



We hope that the description of unconstrained initial data for gravity presented in this 
short article may be useful for such applications like: 

• construction of initial data for numerical analysis (see e.g. |13j). 

• description of the so-called dynamical horizons (see e.g. 

• description of initial data which are sufficiently close to spherical symmetry (one 
can use analysis based on different /3-foliations, analogous to the one presented 
in jS] for (5 = foliation), 

• construction of trapped surfaces (see e.g. 8,), 

and in other cases which are not yet discovered or not known to the authors. 
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